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Abstract 

We perform Hartree-Fock calculations to show that quantum dots (i.e. two 
dimensional systems of up to twenty interacting electrons in an external 

^ ■ 

O ' parabolic potential) undergo a gradual transition to a spin-polarized Wigner 

' crystal with increasing magnetic field strength. The phase diagram and 
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ground state energies have been determined. We tried to improve the ground 



state of the Wigner crystal by introducing a Jastrow ansatz for the wavefunc- 

S" 

tion and performing a variational Monte Carlo calculation. The existence of 



so called magic numbers was also investigated. Finally, we also calculated the 
heat capacity associated with the rotational degree of freedom of deformed 



^ ■ many-body states. 

I. INTRODUCTION 

Quantum dots have been subject of recent intense experimental and theoretical research. 
The interest in those nanostructures arises not only from possible new technological appli- 
cations, but also from the desire to understand the fundamental physical problem of a few 
(m < 300) interacting electrons in an external potential and a strong magnetic field. 

In a weak magnetic field electrons form a rotationally symmetric state by occupying 
the lowest Fock-Darwin levels. With increasing magnetic field, the behaviour is determined 
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by three basic mechanisms: the spin alignment of the electrons (eventually resulting in a 
spin polarized system), the Pauli principle (which causes the ground states of the system to 
occur with certain "magic numbers" of the total angular momentum |I||§), and finally the 
Coulomb interaction. This latter makes the electron droplet susceptible to edge excitations 
and bulk instabilities, as the electron-electron interaction favors a larger area HQ. 

The scenario above assumes a unbroken rotational symmetry. Questions concerning a 
high-field transition to so-called Wigner molecules and crystals, in which electrons occupy 
fixed sites in a rotating frame and are therefore localized, have been investigated by Maksym 
|| and Bolton et. al. ||. Maksym considers a 'large angular momentum limit' of systems 
of up to five electrons and describes only excited states of integer angular momentum. He 
speculates on the existence of ground state Wigner molecules in the large field limit. Bolton 
et. al. simulated up to forty classical interacting point charges in an external parabolic 
potential. 

In this article we consider the ground state properties of up to twenty electrons in the 
limit of a strong magnetic field. We treat the full quantal problem by solving the Hartree- 
Fock equations for the two-dimensional electron gas. After a brief description of the well 
known phases of the rotationally symmetric they appear in the Hartree-Fock ap- 

proximation, we present the gradual transition towards a Wigner molecule and crystal; the 
various spatial configurations are shown and the shell structure is compared to that of the 
classical calculation ||. We further describe a variational Monte Carlo calculation which 
we performed to improve the wave function by a Jastrow ansatz. Finally, we investigate the 
rotational spectra associated with the breaking of the continuous rotational symmetry; the 
heat capacity associated with this new rotational degree of freedom is calculated. 

II. THEORY 

We consider m electrons of effective mass m* in a plane (x, y) confined by an external 
parabolic potential, V(r) = ^m*UQr 2 , and subject to a strong magnetic field B = B e z . The 
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Hamiltonian for such a system is 

m -i mi m, * R C ^2 

* = E ^A 2 + E + rf) + E ^ Li ^ + E (i) 

i=l z/7t j=i z i=i i<j e l'i ' jl 

where ITj = 4Vj + ^(ri) is the kinetic momentum of the ith electron in the vector potential 
A(fl) = ^f(—yi,Xi,0). We include the spin degree of freedom of the electrons, Si (which 
implies the addition of the Zeeman energy with g* pa 0.54 as the effective g factor for the 
material GaAs), but neglect any spin-orbit interaction. (An order-of-magnitude estimate for 
the magnetic field strength induced by the circular motion of an electron with an angular 
momentum of h at a distance of l ~ 1.0 x 10 cm is only — * — ft ~ 2.8xl0~ 5 T.) Defining the 
frequencies u c = and uj(B ) = \Juj 2 + we rewrite the coordinates as dimensionless 



complex variables, Zi = ^ ? where l = y m *^ Bo ) is the magnetic length, to obtain 
m h 2 ( 1 \ h 

H = J2 -A T^A + m*LU 2 (B )l%\ Zl \ 2 + -ujfrd* - Zi d Zi ) 

i=i m Vo ) 1 

H is the single-particle Hamiltonian whose eigenfunctions will form the basis states for 
the Hartree-Fock calculation. For the spatial part of H we define 

a f = —(z-d z ), = -j=(z-d g ), a = -j=(z + d s ), b = -^={z + d z ) (3) 

with [a, a"!"] = [b, tf] = 1 and write 

I i 

H (z) = hu(B )(a!a + tfb + 1) - -hu c {b ] b - a f a) = huj(B )(2a^a + C + 1) - -huj c C, (4) 

where £ := tfb — a^a is the angular momentum of the particle. The eigenvalues of the single 
particle Schrodinger equation H ^/ nk = e n k^ n k, £nk = hu(B ){2n+ \k\ + l} — jhu c k, indicate 
that for strong magnetic fields, e.g. B pa 10 T, one has hu c ~ 17 meV and hu(B ) pa 9 meV, 
and all particles occupy the Fock-Darwin states with n = 0: 

Ae = e n=0 L - e n =0L-i ~ 0.5 meV < e n=0 L - e n =iL ~ 18 meV. 
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We therefore restrict our calculation to the n = level, which resembles the lowest 
Landau level if the external potential had been switched off. The eigenfunctions are 

~ z k exp(—\z\ 2 ) and are identical to the usual form 0^ ~ e~ lfcv r fe Ln (r) exp(— ^), l)n\r) 
being the Laguerre polynomial of degree zero. 

In the Hartree-Fock calculation, we minimize the Hartree-Fock energy 

E HF = ($\H\$) = ^t llh (®\c\ 1 c l M + \ E \hhiA®\4A c h c iM ( 5 ) 

= ^2 t hhPhh + 2 E Phh^hhhhPhh 
hh hhhh 

with p^i 2 = (<&\c\ c^l®) being the density matrix and |<3?) being a Slater determinant, c] 
creates a fermion in the state ^i, while its Hermitian conjugate q destroys it. The indices 
U = (ki, Si) run over all orbital states k, as well as the spin degree of freedom s = {+§, — |}. 
ti x i 2 = \Ho\^i 2 ) = ek 1 6k 1 k 2 + (~ l) Sl+ ^ 9 ^^ Ba 5 SlS2 is the single particle matrix element of 
the Hamiltonian and v the antisymmetrized Coulomb matrix element, 

- _ e 2 I (ablj^lcd) - (ablj^ldc) if s c = s d , 

Vabcd — r- \ \P) 

To minimize (|5]), we vary with respect to p, 5E HF /5p = 0, with the contraints that we 
stay within the set of Slater determinants (p 2 = p) and conserve the number of particles 
(trp = m) resulting in the matrix diagonalization problem 



X hijDjk — E (tij + r»j) Djk — E I tij + X/ Vil'jlPW J Djk — ZkDik, 



(7) 



where V is the so-called mean field. The eigenvectors of h represent the new single- 
particle states {k}, that are to be occupied according to the energies Equation (^) has 
to be solved self consistently, since pu> = YJlLi DuDli- 

We use the Fock-Darwin representation in our calculation and take into account up to 
200 single-particle states (including spin). We tested our code by comparison with the 
results of Pfannkuche et. al. |]] and Bolton [|7|]. Although Pfannkuche et al. describe 
quantum dot helium in a model space that is different from ours (they included n ^ 



states in their calculation), our ground state energies of total angular momentum J = 1 for 
T < B < 5 T coincide with their Hartree-Fock calculation within less then 2%, and the 
J = ground state energies agree with less than 5%. As one can see in table I of reference 
H, the n/0 coefficients in their J = ground state are larger than in their J = 1 ground 
state, so that the n^O space is more significant for those magnetic field strengths. Similar 
results are obtained if we compare our results to the fixed node Monte Carlo calculation of 
reference [0]. In the spin polarized case our ground state energies agree within a few percent, 
while we overestimate the energy of the depolarized system by up to 15%. This is due to 
the larger correlation energies (ignored in a Hartree-Fock calculation) when two electrons 
can occupy the same orbital. Since the questions addressed in this article concern the spin 
polarized regime, this deviation from the results of is of little concern. 

The Hartree-Fock approximation is known to conserve symmetries present in the initial 
trial wavefunction. To generate deformed solutions, we started with a quite arbitrary, but 
not rotationally invariant, initial Slater determinant, which produces a deformed initial 
mean field. Self consistent iteration of the Hartree-Fock scheme guarantees amplification 
of solutions with the symmetry of the Wigner molecule. Of course, the same converged 
solution must be reached for several different initial states to give confidence that it is the 
true minimum. 

III. NUMERICAL RESULTS 

We have used the material constants of GaAs (i.e., m* = 0.067 m e and e = 12.9, as well as 
an external potential strength of huj = 3 meV) for our calculation. To observe the expected 
phase transition, we first consider a system of m = 10 electrons. Fig [1] shows the ground 
state energy as a function of the magnetic field strength B and for comparison the lowest 
energy of the rotationally symmetric system. The Wigner molecule becomes the ground 
state for Bq > 5.2 T, while at smaller strengths the rotationally symmetric state is favored. 
The system undergoes spin polarization from B = T to 1.5 T, where the spin polarized 
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so-called maximum density droplet || prevails. At B$ = 4.5 T bulk instabilities result in 
unoccupied inner Fock-Darwin states. The transition to a Wigner molecule, and later to 
a crystal, happens very gradually. (We refer to the case where the probability density is 
deformed, but still very smeared out as a "molecule" , while a "crystal" signifies well localized 
and distinguishable electrons, as illustrated in Fig ^|.) The molecule at B = 6 T is lower 
in energy by only 0.2% (0.542 meV) relative to the rotationally symmetric solution, while 
the crystal at Bq = 10 T gains about ~ 3 meV, which is of order of the strength of the 
confining potential. Note that the deformed ground states are not eigenstates of the total 
angular momentum operator J = 

The rotationally symmetric case suffers further complication with increasing magnetic 
field strength: While at first (B « 6 T) the hole in the bulk widens (the / = 1, 2, 3 
Fock-Darwin levels empty) and later at Bo = 6.75 T a fourth state empties, the solution 
transforms into two separate rings at Bo > 9 T. 

For further insight into the various transitions, the radial (angle- averaged) particle dis- 
tribution for various magnetic fields is shown in Fig [5] for m = 20 electrons. The crystalline 
state has one electron in the center of the dot, seven in a middle ring and 12 electrons in the 
outmost region. Correspondingly, the B = 20 T curve shows three maxima. For B = 6 T 
the center electron and the seven in the middle ring have almost uniformly merged to a flat 
distribution which extends to z ~ 2, and the outer ring can now be found at z ~ 3. For 
B = 3 — 4 T we find again the so-called maximum density droplet: the electrons occupy the 
first twenty Fock-Darwin levels, since they are polarized. Further lowering of B results in 
a depolarization, allowing further accumulation of electrons near the origin. Since we only 
take into account n = states, we cannot claim to represent the physical situation for the 
smaller field strength, although we do reproduce the energies in this regime quite well, as 
noted above. 

In Fig ^ we plot the separation energy A(m) = E m+ i — E m and the differences in the 
separation energy, A 2 (m) = A(m + 1) — A(m), as functions of the particle number m in the 
crystal regime, Bq = 20 T. There is a large drop in A 2 of ~ 0.5 meV whenever charge can be 
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put to the outer region of the dot (see, e. g., m = 4 and m = 8), in accord with charge being 
distributed over a larger area, thereby reducing the Coulomb energy. In the case where one 
charge is placed in the center and two rings outside (m = 14), the gain in energy is reduced 
by the fact that more particles outside feel a stronger external potential. The tendency here 
is that the Coulomb energy plays a less and less important role, weakening the slope in 
the separation energy, combined with the fact that one can pack more particle in the outer 
region. 

For comparison, we also show in Fig || A(m) and A2(m) for the lowest rotationally 
symmetric state. No clear tendency in the behaviour of A2(m) is evident. The system is 
frustrated by the particles having to occupy Fock-Darwin levels. 

In Table | we show the spatial configurations of the system in the Wigner-like structure 
(obtained by enumerating the number of electrons occupying the corresponding rings) and 
give the ground state energies. We generally confirm the configurations of the classical 
calculation of reference || as well as the exceptional behaviour of the m = 6, 10, 12 and 
17 clusters, although there is no peak in A 2 for m = 14 (the peaks in Fig [| correspond to 
the cusps of Fig 5 of reference 0), since ten electrons are moved outside for the m — 16 
configuration. 

The Hartree Fock calculation is based on a theory of independent particles moving in an 
average potential. We improved the wavefunctions for the Wigner regime to a many body 
wavefunction by introducing a Jastrow type function 



\ l / 

where S is the symmetrizer and &hf the Hartree Fock solution to the problem. To guarantee 
a convenient symmetrized form of the product of these function, we made the ansatz 



for the pair correlation function f(z{ — Zj) with k as a variational parameter. We performed 
a variational Monte Carlo calculation || to evaluate the energy 
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(8) 



fyZi Zj^j \Z{ Zj 



(9) 



E \h\ — /[in ■ (10) 

As it turned out, the Jastrow type wavefunction did not significantly improve the Hartree- 
Fock energy. In the case of 10 electrons and B = 20 T, the energy could only be improved 
by 0.1% (0.4 meV) at k = 0.1. For 0.1 < k < 1 the energy is slowly increasing, while for 
k > 1 highly excited states are simulated as more holes are introduced into the wavefunction. 
Obviously, the Hartree-Fock solution already describes the Wigner state accurately. 

In Fig we plot the phase diagram with respect to number of particles and the ratio 
We omit the regime of bulk instabilities, since it is of minor importance. The phase 
boundary of the spin polarized regime and the partially unpolarized regime suffers again 
from the Hartree-Fock approximation, as it bends down with decreasing number of electrons. 
The boundary of the molecular regime is defined by how much the continuous rotational 
symmetry is broken: the fractional uncertainty in the total angular momentum is / = jjjr = 
A ^ — , and we define a molecule by / > 1%. The boundary is fairly constant for m > 6, 
but, since the transition is gradual, it has some uncertainty. For less than eight particles, 
we find a small drop in the boundary, due either to some non-obvious physical effect or to 
the approximation we use. 

In our Hartree-Fock solutions of ten or more electrons and B = 20 T, the relative 
uncertainty in total angular momentum, /, is of order of 10%. As in atomic nuclei, these 
deformed solutions give rise to rotational spectra, which do not appear in the case of the 
unbroken symmetry. We have estimated the spectrum of rotational excitations by projecting 
the Hartree-Fock Slater determinant onto eigenfunctions of good angular momentum / || . 
The projector has the form 

pj = i r 2 « e ia(j-i) da 



2tt Jo 

and the energies which result from taking the mean value of H with the projected wave- 
functions are given by 



^ _ (^HP 1 ^) _ Jdah(a)e- ihl 
($|P J |$) jdan(a)e- Uo 



^proj - ,^,An^ ~ fJ„.„/„.L->Tr,' [l -> 



s 



defining the quantities h(a) = (<&\He iaJ \§) and n(a) = ($|e iaJ |$). Since the standard 
deviation in J is only of few percent, one can calculate these matrix elements approximately 
by writing h(a) in the expansion 



h(a) 



+ ^) n(a, (13) 



One justifies this ansatz with the fact that it represents a Taylor expansion of the Fourier 
transformed function h(a)/n(a), and, assuming that both quantities are sharply peaked at 
a = 0, this quotient is smooth and can be approximated by a few terms of equation (|13|). 
By operating ((—J) + f gs) on equation ( |T3|) and setting a = 0, one gets an inhomogeneous 
system of equations for the unknown ho...h,K' 



A 



(tf (AJ) m > = 5> n ((AJ) m+n >. (14) 



n=0 



Equation (|T2") can then be expressed as 



A" 



trE^-Wf. (15) 



n=0 



We restrict ourselves to K = 2, since higher terms involve the calculation of fc-body operators 
with k > 4. For this case, we have 

{AJ 2 )(HAJ 2 } - (AJ 3 ){HAJ) - {AJ 2 ) 2 (H) 
2 ~ (AJ 2 )(AJ 4 ) - (A J 3 ) 2 - (A J 2 ) 3 

(AJ 2 ) (AJ 2 ) 
h =(H)-h 2 (AJ 2 ). (16) 

Fig [| shows the rotational spectra for m = 10 and m = 20 electrons and B = 20 T as 
a function of the quantum number /, where we have substracted the shifted ground state 
energy, which is obtained from the Hartree-Fock energy (H) by substracting the spurious 
rotational energy /i2(AJ 2 ), which is only of order 0.25 meV in both cases. The moments 
of inertia associated with these states are 

J Y = JL = 5.2 x 10 5 h 2 /eV for m = 10 and 

J y = 1.9 x 10 6 h 2 /eV for m = 20. 

9 



In order to excite a molecule with circular polarized radiation, one has to produce photons 
of minimal energy of AE^^I = 224) = E* roj (I = 224) - ((H) - h 2 (AJ 2 )) = 1.12 • 1(T 7 eV 
for the ten electron molecule and AE^ 20 '(I = 790) = 3.2 ■ 10~ 8 eV for twenty electrons, 
which are the energy differences between ground and first excited state. These energies 
correspond to radio frequencies of = 27.06 MHz and z/ 20 ) = 7.73 MHz. Note that the 
corresponding wavelengths are in the transparent region for GaAs. Therefore the measure- 
ment of transmission coefficients of circular polarized radiation should give experimental 
evidence of Wigner molecules. The level spacing, AE « — t^-AI, of the excited states 
then increases with higher states, resulting in excitations in the microwave region. The heat 
capacity connected with this rotational degree of freedom, 

d<U> 3 lL„, / EL n ; 



c 



E^exp , (17) 



dT dT Z VT W ° 3 V k sT 



where Z = l + X)/exp (^—Ep ro j/kBT^j is the partition function and ks Boltzmann's constant, 
should therefore reach its classical value of even for temperature as low as 1 K. Fig [7] 
shows the well known Schottky anomaly of the heat capacity, typical for a system where only 
two states are of importance, at low temperatures of ~ 1 mK. As expected, it approaches 
for high temperatures. For the indicated temperature regime the heat capacity has 
converged within our model space, which consists of 400 rotational states and shows the 
expected typical behaviour of a quantum mechanical rotor in a heatbath. 

The energy levels of the vibrational modes of a single electron in the crystal can be 
estimated in a simplified one dimensional model. Concerned only with the radial degree of 
freedom, an outer electron (in the case of ten electrons) interacts with the external potential 
and the Coulomb potential of the two inner electrons, which we regard positioned at the 
center: 

1 2e 2 

V(r) = -rrfulr 2 + . (18) 

Expanding the potential around the equilibrium position r$ of the outer electron to second 
order, we obtain 
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../ , 1 „. f o 86 I . >2 f * 2 26 I . . 1,22 ^6 / ir ,\ 

\/(r) = -m ^ + — — Mr - r ) + 4 m u r ^ r-r +-m ^ r + . (19) 

2 [ em*rg J [ erg J 2 er 

The electron is confined by the parabolic part of this expansion with an corrected strength 
uj' = ^JuJq + ^f^-- Setting r ~ 2 x 1CT 6 cm, the energy levels for the vibrational modes of 
the electron are separated by AE vib = hu(B ) = ^ (hu 1 ) 2 + \(huj c ) 2 « 21 meV, much larger 
than the separation in the rotational energy levels (AE rot w 10~ 4 meV around / = 225 h). 
Vibrational modes therefore contribute only marginally to the heat capacity and can be 
easily suppressed by proper excitation of the rotational modes only. 

In summary, we have shown in a full quantum mechanical treatment that there exist 
regimes where Wigner molecules and crystals are the ground states of quantum dots. We 
have also described rotational spectra of quantum dots, which arise from the existence of 
deformed Hartree-Fock solutions. This broken symmetry could make it possible to detect 
Wigner molecules experimentally by exciting the rotational excited states of the system. 
Open questions remain: How do the two-body correlations neglected in the Hartree-Fock 
approximation influence systems of few electrons (m < 5) and what kind of state is formed 
in the case of many particles (m > 40)? 
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TABLES 



number of 
electrons 


energy 
[meV] 


ring occupations 
inner — middle — outer 


1 


17.247 


1-0-0 


2 


40.085 


2-0-0 


3 


66.439 


3-0-0 


4 


96.463 


4-0-0 


5 


129.986 


5-0-0 


6 


166.346 


1-5-0 


7 


205.448 


1-6-0 


8 


247.636 


1-7-0 


9 


292.621 


2-7-0 


10 


339.934 


2-8-0 


11 


389.489 


3-8-0 


12 


441.634 


3-9-0 


13 


496.008 


4-9-0 


14 


552.825 


4-10-0 


15 


611.879 


5-10-0 


16 


673.004 


1-5-10 


17 


736.135 


1-5-11 


18 


801.162 


1-6-11 


19 


868.558 


1-6-12 


20 


937.973 


1-7-12 



TABLE I. Ground state energies and spatial distributions of Wigner crystal in quantum dots 
for up to twenty electrons at B = 20 T. 
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FIGURES 

FIG. 1. Ground state energy of the Wigner molecule (solid line) and lowest rotationally sym- 
metric state (dashed line) for ten electrons. 

FIG. 2. The electron density distribution for m = 10 electrons. The left panel shows the 
solution at Bo = 6 T, the right at Bq = 10 T. 

FIG. 3. Radial electron density of m = 20 electrons for different values of Bq. The solid curves 
represent Wigner crystals, while the dashed curves show the slow transition to a maximum density 
droplet, which is drawn with a dashed-dotted line. Depolarization sets in for the cases of the dotted 
curves. 

FIG. 4. Separation energy A(m) and the differences in the separation energy A^m) for Wigner 
molecules (upper two diagrams) and for the lowest available rotationally symmetric states (lower 
two diagrams). 

FIG. 5. Phase diagram for quantum dots, plotting ^ versus the number of electrons (m). The 
lines crudely trace the boundaries to guide the eye. 

FIG. 6. Rotational spectra for ten electrons (upper diagram) and twenty electrons (lower di- 
agram), when they have formed Wigner molecules at Bq = 20 T, as a function of total angular 
momentum I 

FIG. 7. Heat capacity c arising from the rotational spectra of Figure ||. The dashed line shows 
the twenty electron system, the solid one the ten electron case. 
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